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TIME-DEPENDENT SLOW-TIME-SCALE THEORY 
OF FREE-RUNNING AND PHASE-LOCKED 
GYROTRON OSCILLATORS 

I. INTRODUCTION 

There is currently considerable interest in the development 
of high power phase-locked gyrotron oscillators. These devices 
have the potential to combine the high efficiency and power 
associated with oscillators with the coherence and phase control 
properties associated with amplifiers. Although previous 
theoretical work on gyrotron operation[1-8] has been successfully 
applied to the development of cw devices for heating of tokamak 
plasmas, the investigation of phase-locked gyrotron operation, 
particularly involving pulse-line accelerator driven devices, 
presents a new challenge to the theorist. These devices, based 
on intense relativistic electron beams, are characterized by 
short pulselengths and voltage risetime and ripple effects which 
can prevent the establishment of a steady-state equilibrium. A 
consideration of non-stationary operation is also involved in the 
investigation of conditions for phase-locked operation. Time- 
dependent effects can be studied using a particle-in-cell 
simulation code of the type developed by Lin and co-workers[9]. 

In this work an alternative slow-time-scale approach is used to 
obtain a relatively simple, economical analysis of the time- 
dependent, nonlinear dynamics of pulsed gyrotrons. 

The time-dependent theory of gyrotrons has been considered 
by Nusinovich and co-workers[10,11], mainly in the context of 
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multimode operation and mode stability. A time multimode theory 
of quasi-optical gyrotrons has been developed by Bondeson et 
al.[12]. Early work on mode selection and phase-locking of 
vacuum tube oscillators was carried out by Van der Pol[13] and 
Adler[14]. An analytical theory of the conditions for phase¬ 
locking gyrotrons has been presented by Manheimer[15,16]. The 
emphasis in the present work is on the modeling of transient 
effects in gyrotrons whose pulse length is only a few times 
longer than the characteristic risetime of cavity fields. The 
approach is developed for the case of several interacting modes 
with closely spaced frequencies and interacting at the same 
harmonic of the cyclotron frequency. Calculations are presented 
for single mode operation of relativistic gyrotron with 
parameters similar to an experimental device under test at the 
Naval Research Laboratory[17]. Both free-running and phase- 
locked operation are considered. The phase-locking signal is 
introduced via a prebunched beam. The effect of a voltage ripple 
is investigated. The possibility of phase-locked operation in 
the presence of voltage ripple - an important issue for proposed 
NRL experiments - is demonstrated. Calculations for multimode 
operation will be the subject of future work. 

The time theory of gyrotron oscillators derived in this work 
is based on slow-time-scale equations for the electron motion 
similar to those used in steady-state models[7,18,19]. Slow- 












time-scale equations for the cavity rf field amplitude are 
obtained by expressing the time-dependant behavior relative to a 
reference frequency w 0 which is close to the operating frequency 
<t). The next section of this paper contains an outline of the 
theoretical approach. The third section contains the results of 
calculations and the last section presents some conclusions drawn 


from this research. 




















II. THEORY 


Consider a gyrotron with a cylindrical resonator and a thin 
annular beam. The electrons follow helical trajectories in the 
applied axial magnetic field about guiding centers located at a 
radius R 0 from the symmetry axis. The electron beam interacts 
with one or more resonator modes, which are assumed to be TE 
modes with closely spaced oscillation frequencies. It is 
convenient to look at time effects which remain after a reference 
frequency w 0 has been factored out. These effects are 
characterized by time scales which are much longer than the wave 
period and are incorporated in time mode amplitudes A n (t). The 
total transverse electric field is expressed in the form 


E t - Z A n (t )h n ( z )e n ( r , 9; z )exp<-ioo 0 t) (1) 

n 3 1 

where e n is a transverse mode vector function defined in Ref. 18, 
and h n (z) is an axial profile function which satisfies the 
Helmholtz equation; 










The function P n »h n e n is an eigenfunction of the empty resonator 
in the single transverse mode approximation with eigenfrequency 
w no . Thus the present formulation neglects beam loading effects 
on the mode axial profile. N is the number of interacting modes 
The slowly varying mode amplitudes have the form: 

A n (t) - f„ me -1 *" (t) (3) 

The phase derivative ^J-d^/dt corresponds to the frequency 
difference between the mode oscillation frequency « n and « 0 . 

The transverse electric field satisfies the wave equation: 

. 3 2 E t 3 J t 

72 Et _ A- =. u 0 - (4) 

c 2 3t 2 3t 

where J t is the transverse ac current density, c is the speed of 

light and fj a is the permeability of free space. MRS units are 
used throughout except as noted. Substituting Eq.(l) into 
Eq.(4), introducing a mode Q factor, Q n , to account for 
diffraction and ohmic losses, and noting that dA n /dt<<w 0 A n one 


obtains 


















where 40 no is the empty resonator eigenfrequency for the nth mode. 
Multiplying Eq.(5) on both sides by (h n e n )*, integrating over the 
cavity volume and applying the orthogonality properties of the 
transverse vector functionsl18] leads to the mode excitation 
equation: 


r d 

W o 

- i c 2 fj Q - 


+-i(w 0 

~ “n.lkltl " dzda 

L dt 

2Q„ 

J 2N n c o i 

Nn 

is the axial 

mode normalization integral: 


iu> t 


N n - Pdz |h n (z) 

J o 


To evaluate the right hand side of Eq.(6) the rf current density 
is written in the approximate forra(7]: 


J,. =J exp ( -ioo t) 
4 


where 


J - f d(u 0 t)J t e 1 
w o J n 


iw 0 t 


Introducing the normalized frequency shift 6 n = ( -oo n o )/t»> 0 and 
time t=<a) t and substituting Eq. ( 8 ) into Eq. ( 6 ) leads to: 


r d 1 i c _ 

— +-i6 n A n (x) -- dzda h; e; • J e - 1 

L dx 2Q J 2N_w . J 




























To calculate the ac current density, the interaction with the 
electron beam is treated in the single particle approximation. A 
considerable simplification of the general time-dependent problem 
results if one uses the fact that the characteristic rise-time of 
fields in the resonator is much longer than the electron transit 
time in the cavity. In this case one can use a quasi-steady- 
state approximation in which the electron trajectories are 
calculated for rf fields with fixed amplitude, : and 

linearized phase, '(t 0 )(t-t 0 ). The slow-tiroe-scale 

nonlinear electron equations of motion for a particular harmonic 
are readily deduced from previous steady-state analyses[7,18] and 
are given by: 
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where u t «yv t /c is the normalized transverse momentum amplitude, 
u z -rv,/c is the normalized axial momentum, A is a slowly varying 
transverse momentum phase, s is the harmonic number, m n is the 
mode azimuthal index, y is the relativistic mass ratio, k n is the 
mode transverse wavenumber, r L is the Larmor radius of the orbit, 
J s (J' s ) is (the derivative of) a regular Bessel function, 5 is 
the relativistic cyclotron frequency correspondong to the input 
beam energy, E 0 is the orbit guiding center azimuthal angle, and 
F n is the normalized rf field amplitude: 


F n " 


J (k n Rj 
m n -s n ° 


m.c 2 /[ n( l-m2/x^ 2 ) ] J ( x^ ) 

m n 


( 12 ) 


Quantities with a """ have been normalized according to: 

£-z/r uo , ? L -r L /r wo , C-Cr M0 /c, ^ 0 -« o r wo /c, and & n -k n r uo . R o 

denotes the orbit guiding center radius, e is the electron 
charge, m # is the electron mass, x' n is a zero of J' for the mode 
n, and r wo is an arbitrary normalization length usually chosen to 
be a representative cavity wall radius. 

The ac current density is obtained by integrating Eqs.(ll) 
for an appropriate set of initial conditions at the cavity input 
at z 0 . For a cold, phase-mixed electron beam: u t (z 0 )=u to , 
u,(z 0 )-u #0 , and A(z 0 )=A 0 is uniformly distributed in the interval 
[0,2n]. The case of an electron beam prebunched by passage 
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through a bunching cavity driven at a frequency « Q can be modeled 
by an initial phase angle distribution of the form: 


A 0 - 0 o + q B sin0 o - * 0 (13) 

where 0 O is uniformly distributed in the interval [0,2rt], q B is 
the bunching parameter, and i 0 is an initial phase. For an 
annular beam, Eqs.(ll) are averaged with respect to S 0 . For a 
thin annular beam the transverse ac current density is given by 



Substitution of Eqs.(14) and (9) into Eq.(10) leads to the 
following equations for the mode amplitude and phase: 


dr 


Re 
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* - - i ( A+'l'_ - ( m -s ) E 


dzh; <— J' ( k n r L )e 



(15a) 


d'J'n 

dx 


6 n - im 
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dzh; <— J' ( k n r L )e l[ A+ * n ( S) “° ] > A - 


v. 

u z 0 ' “o . 


(15b) 


where I is the normalized current parameter: 
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r wo J m-s (k n R o } 


m 0 ca> 0 n(l-m 2 /x' 2 ) J 2 (x' )N r 


(16) 


and < > a denotes an average with respect to the variable a. The 
time behavior of a free-running or driven gyrotron oscillator is 
obtained by solving Eqs.(lS) together with (11) subject to 
appropriate initial conditions as discussed above. When the 
phase \|/ n of a driven oscillator (via beam prebunching) converges 
to a constant value, the oscillator is said to be phase-locked. 
The locking frequency is &> 0 in the present formulation. When the 
interaction occurs with only a single circularly polarized mode, 
the ac current density averaged with respect to A 0 is independent 
of the guiding center azimuthal position so that the average with 
to S 0 may be omitted. 
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III. RESULTS 


Calculations have been carried for a high voltage 35 GHz 
gyrotron configuration similar to the NRL experiment recently 
reported by Gold et al.(17]. The peak voltage and current of the 
annular beam are taken to be 650 kv and 1.5 kA, respectively. 

The beam guiding center radius is R 0 -1.16 cm, the cavity radius 
is 1.6 cm, and the longitudinal profile of the cavity fields is 
assumed to be sinusoidal with length 4 cm. The operating mode is 
the TE g 21 circular mode with polarization counter-rotating to the 
beam rotation. The beam pitch ratio asv t /v 2 =>l. The cavity Q 
factor Q«250. The applied magnetic field is 2.4 Tesla. Spreads 
in beam guiding center and pitch ratio are neglected as are 
space-charge effects. In these calculations the cavity is seeded 
with low amplitude rf fields (1 kV/cm at the beam in the present 
case). The results of calculations for free-running and phase- 
locked oscillator operation are presented below. The effect of a 
voltage ripple is shown. 

The first set of calculations correspond to an idealized 50 
nsec voltage pulse with a 3 nsec rise from 75% of peak to peak 
voltage and no ripple. The beam current and pitch ratio are 
assumed to vary linearly with voltage and as V 1 - 33 , respectively, 
during the voltage rise up to the maximum values given above. 
Results for the free-running oscillator are shown in Figures 1 
(a)—1(c). These Figures show, respectively, the electronic 





















efficiency, operating frequency, and rf phase of the oscillator 
plotted as a function of time. Figure 1(b) plots the shift in 
operating frequency w relative to oo 0 which is taken to be the the 
cold cavity eigenfrequency (which is close to 35 GHz). These 
plots indicate that steady-state operation is achieved after 
about 8 nsec. The steady-state efficiency is 17.7%, 
corresponding to an output power of 173 MW, and the steady-state 
frequency shift 6 r R «(aj r R -co 0 )/w 0 due to beam loading is 5.5xl0 -4 
which corresponds to 19.2 MHz. This frequency shift leads to a 
linearly increasing rf phase with time as shown in Figure 1(c). 

Results for an oscillator driven via a prebunched beam are 
given in Figures 2(a)-2(c). The beam is modulated at a frequency 
oo 0 and the modulation amplitude is modeled by the bunching 
parameter q B =0.5. An estimate for the locking signal power 
required to obtain a beam with q B »0.5 is given below. In this 
calculation the difference between the normalized locking signal 
frequency and the normalized free-running oscillator (FRO) 
frequency is S L = ( -<*) F R ) /oo 0 =*1.0 5 xlO~ 3 or Af=35 MHz. Figures 2(b) 

and 2(c) show that phase-locked operation is obtained after 30 
nsec into the voltage pulse. The phase-locked condition is 
indicated by a frequency shift go-<*> o =0 and constant rf phase. As 
shown in Figure 2(a), in this case the oscillator efficiency is 


reduced slightly to 16% while phase-locked. 













If the difference between the locking signal frequency and 
the FRO frequency Af is too large, phase-locking will not occur. 
Instead a periodic perturbation of the oscillator behavior is 


observed as shown in Figures 3(a)-3(c). The locking signal 
detuning in this case is S L »1.55xl0 -3 or Af-54 MHz. The beat 
frequency of the perturbation is about 28 MHz. The non-phase- 
locked behavior of the rf phase is clearly shown by Fig. 3(c). 

The case of a FRO with a voltage ripple is shown in Figures 
4(a)-4(c). The voltage ripple amplitude is 2% and the ripple 
period is 10 nsec. Other parameters are the same as for the FRO 
case discussed above. Figures 4(a) and 4(b) show that this 
voltage ripple causes a 3.4% amplitude ripple in the efficiency 
and a 0.06% amplitude ripple (21 MHz) in the frequency. Figures 
5(a)-5(b) show the results of driving this oscillator 
configuration with a prebunched beam with q„=0.5 and modulation 
frequency shift & L =5.5xlO -4 (19 MHz). Figures 5(b) and 5(c) show 

that phase-locked operation is achieved in the sense that after 
about 30 nsec the oscillation frequency averaged over the ripple 
period equals the drive frequency and the averaged rf phase angle 
is a constant. As in the FRO, the voltage ripple leads to a 
0.06% frequency ripple in the PLO. The frequency ripple causes a 
corresponding 13° amplitude ripple in the phase angle. Figures 
6(a)-6(c) represent the case of an unphase-locked, driven 
oscillator with a 2% voltage ripple. The drive signal frequency 
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shift is S l -2.5x 10“ 3 (88 MHz). The beat frequency for the 
perturbed oscillator is about 34 MHz. The failure to achieve 
phase-locked operation is clearly shown by the phase angle 
plotted in Figure 6(c). 

To estimate the input power needed to obtain a beam bunching 
parameter q B -0.5 for circularly polarized radiation and a 
gaussian rf profile in the bunching cavity, one can apply the 
small signal gyroklystron theory developed by Tran et al (20]. 
The required input power in Watts is: 


P in (W) 


^c 3 m 2 (x^ 2 -m 2 ) J 2 (x' n ) y| 6? 0 3 Z0 q| e‘ 
2’I 2// 0 e 2 J 2 (k n R 0 ) (4 3 Q 1 


(17) 


where /y x -nH 2 to L 1 /(3 zo X is the normalized interaction length of the 
bunching cavity, u A is the normalized drift length, Q x is the 
bunching cavity Q factor, and x= ( w Q -S2) L x /4 v 2 o the bunching cavity 
magnetic field detuning parameter. The other parameters are 
defined above. Typical choices for ij 1 and fj A are 2 and 3, 
respectively. Assuming the same mode, magnetic field and beam 
parameters in the bunching and power cavities, and Q 1 =200, one 
obtains P in =*530 kw as the power required to phase-lock a gyrotron 
oscillator with an output power of over 160 MW for an effective 
gain of over 24 dB. 
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IV. CONCLUSIONS 


A time slow-time-scale theory and a computer code for 
numerical calculations have been developed to investigate pulsed 
free-running and phase-locked gyrotrons. This code is expected 
to be a valuable resource for the design and interpretation of 
planned phase-locked gyrotron experiments. The calculations 
presented here show the feasibility of observing phase-locked 
operation with 50 nsec pulses. Phase-locking over a significant 
bandwidth is shown to be possible with a locking power 24 dB less 
than the oscillator output power for a single prebunching cavity. 
The effect of voltage ripple has been investigated and shown not 
to prevent phase-locked operation, however, applications 
requirements for rf phase control may place significant limits on 
allowable voltage ripple. 






















V. ACKNOWLEDGEMENTS 


Helpful discussions with Drs. Steven Gold and John Burke are 
gratefully acknowledged as is the computational assistance of 
Lisa Lagrand. This work was supported by the Office of 
Innovative Science and Technology of the Strategic Defense 
Initiative Organization under a contract administered by the 
Harry Diamond Laboratories. 


» *• A “» -V -V •»».> ■ » - - »« /• A AA A A AAlAJVA.% A AA AlN- 















REFERENCES 


a) JAYCOR, Inc., Vienna, VA 22180. 

b) Electrical Engineering Department, University of Maryland, 
College Park, MD 20742. 

1. P. Sprangle and W.M. Manheimer, Phys. Fluids lj), 224 (1975). 

2. E. Ott and W.M. Manheimer, IEEE Trans. Plasma Sci. PS-13 , 1 
(1975) . 

3. V.A. Flyagin, A.V. Gaponov, A.V. Petelin, and V.K.Yulpatov, 
IEEE Trans. Microwave Theory Tech. MTT-25 , 514 (1977). 

4. P. Sprangle and A.T. Drobot, IEEE Trans. Microwave Theory 
Tech. MTT-25 , 528 (1978). 

5. K.R. Chu, Phys. Fluids 2_1, 2354 (1978). 

6. K.R. Chu, A.T. Drobot, H.H. Szu, and P. Sprangle, IEEE Trans 
Microwave Theory Tech. MTT-28 , 313 (1980). 

7. A.W. Fliflet, M.E. Read, K.R. Chu, and R. Seeley, Int. J. 
Electron. 5_3, 505 (1982). 

8. B.G. Danly and R.J. Temkin, Phys. Fluids .29, 561 (1986). 

9. M. Caplan, A.T. Lin, and K.R. Chu, Int. J. Electron. 5^, 659 
(1982) . 

10. M.A. Moiseev and G.S. Nusinovich, Isvestia VUZ . Radiofizika 
17 , 1709 (1974) [Radio Phys. Quant. Electron. 17_, 1305 

(1976) ] . 


17 





















11. N.S. Ginzburg, G.S. Nusinovich, and N.A. Zavolsky, Int. J. 
Electron. 61_, 881 (1986). 

12. A. Bondeson, W.M. Manheimer, and E. Ott, Phys. FLuids 26, 285 
(1983) . 

13. B. van der Pol, Phil. Mag. 3, 65 (1927). 

14. R. Adler, Proc. Inst. Radio Engrs. 3_4' 351 ( 1946 ). 

15. M.E. Read, R. Seeley, and W.M. Manheimer, IEEE Trans. Plasma 
Sci. PS-13 , 398 (1985 ) . 

16. W.M. Manheimer, Int. J. Electronics 6_3, 29 (1987). 

17. S.H. Gold, A . W . Fliflet, W.M. Manheimer, R.B. McCowan, W.M. 
Black, R.C. Lee, V.L. Granatstein, A.K. Kinkead, D.L. 
Hardesty, and M. Sucy, Phys. Fluids (to be published). 

18. A.W. Fliflet, Int. J. Electron. 61, 1049 (1986). 

19. T.M. Antonsen, W.M. Manheimer, and B. Levush, Int. J. 
Electronics 61., 823 ( 1986). 

20. T.M. Tran, B.G. Danly, K .E.Kreischer, J.B. Schutkeker, and 
R.J. Temkin, Phys. Fluids 29, 1274 (1986). 




































































































































































4740 DISTRIBUTION LIST 
(for 4740 memos) 


Air Force Avionics Laboratory 
AFWAL/AADM-1 

Wright/Patterson AFB, Ohio 45433 
Attn: Walter Friez 


Air Force Office of 
Scientific Research 
Bolling AFB 

Washington, D.C. 20332 
Attn: H. Schlossberg 


Air Force Weapons Lab 
Kirkland AFB 

Albuquerque, Mew Mexico 87117 
Attn: Dr. William Baker 


Columbia University 

520 West 120th Street 

Department of Electrical Engineering 

New York, N.Y. 10027 

Attn: Dr. S.P. Schlesinger 


Columbia University 
520 West 120th Street 
Department of Applied Physics 
and Nuclear Engineering 
New York, New York 10027 
Attn: T.C. Marshall 


Cornell University 

School of Applied and Engineering Physics 
Ithica, New York 14853 
Attn: Prof. Hans H. Fleischraann 
John Nation 
R. N. Sudan 


Dartmouth College 
18 Wilder, Box 6127 
Hanover, New Hampshire 03755 
Attn: Dr. John E. Walsh 


Department of Energy 
Washington, D.C. 20545 
Attn: C. Finfgeld/ER-542, GTN 
T.V. George/ER-531, GTN 
D. Crandall/ER-55, GTN 


1 copy 


1 copy 


1 copy 


1 copy 


1 copy 


1 copy 
1 copy 
1 copy 


1 copy 


1 copy 
1 copy 
1 copy 






MV 

1 









.>v 







■V.li'.l. A .H .!■ A .'I rfl .■« ■■t VL/tyjJlj tj’ki'U Ij'Ll lJ l.li.l U M tJ.MAI.'.; 


Defense Advanced Research Project Agency/DEO 
1H00 Wilson Blvd. 

Arlington, Virginia 22209 
Attn: Dr. S. Shey 

Dr. L. Buchanan 


Defense Communications Agency 
Washington, D.C. 20305 
Attn: Dr. Pravin C. Jain 

Assistant for Communications 
Technology 


Defense Nuclear Agency 
Washington, D.C. 20305 
Attn: Mr. J. Farber 
Maj. J. Benson 
Capt. D. Stone 
Mr. Lloyd Stossell 


Defense Technical Information Center 
Cameron Station 
5010 Duke Street 
Alexandria, Virginia 223m 


Georgia Tech. EES-E0D 
Baker Building 
Atlanta, Georgia 39332 
Attn: Dr. James J. Gallagher 


Hanscomb Air Force Base 
Stop 21, Massachusetts 01731 
Attn: Lt. Rich Nielson/ESD/INK 


Hughes Aircraft Co. 

Electron Dynamics Division 
3100 West Lomita Boulevard 
Torrance, California 90509 
Attn: J. Christiansen 
J.J. Tancredl 


KMS Fusion, Inc. 

39*11 Research Park Dr. 

P.0. Box 1567 

Ann Arbor, Michigan *48106 
Attn: S.B. Segall 


Lawrence Livermore National Laboratory 
P.0. Box 808 

Livermore, California 9*1550 
Attn: Dr. D. Prosnitz 

Dr. T.J. Orzechowski 
Dr. J. Chase 




•\* m w ./ • * 1" O •/ % V 


1 copy 
1 copy 


1 copy 


1 copy 
1 copy 
1 copy 
1 copy 


2 copies 


1 copy 


1 copy 


1 copy 
1 copy 


1 copy 


1 copy 
1 copy 
1 copy 




& 




V-V.A.'.V.V.VV.V.V.V.V.V, 









Los Alamos Scientific Laboratory 
P.0. Box 1663, AT5-827 


Los Alamos, New Mexico 87545 

Attn: Dr. J.C. Goldstein 

1 copy 

Dr. T.J.T. Kwan 

1 copy 

Dr. L. Thode 

1 copy 

Dr. C. Brau 

1 copy 

Dr. R. R. Bartsch 

1 copy 

Massachusetts Institute of Technology 
Department of Physics 

Cambridge, Massachusetts 02139 

Attn: Dr. G. Bekefi/36-213 

1 copy 

Dr. M. Porkolab/NW 36-213 

1 copy 

Dr. R. Davidson/NW 16-206 

1 copy 

Dr. A. Bers/NW 38-260 

1 copy 

Dr. K. Kreischer 

1 copy 

Massachusetts Institute of Technology 

167 Albany St., N.W. 16-200 

Cambridge, Massachusetts 02139 

Attn: Dr. R. Temkin/NW 14-4107 

1 copy 

Spectra Technologies 

2755 Northup Way 

Bellevue, Washington 98004 

Attn: Dr. J.M. Slater 

1 copy 


Mission Research Corporation 
Suite 201 

5503 Cherokee Avenue 
Alexandria, Virginia 22312 
Attn: Dr. M. Bollen 

Dr. Tom Hargreaves 

Mission Research Corporation 
1720 Randolph Road, S.E. 
Albuquerque, New Mexico 87106 
Attn: Dr. Ken Busby 

Mr. Brendan B. Godfrey 

SPAWAR 

Washington, D.C. 20363 
Attn: E. Warden 

Code PDE 106-3113 

Naval Research Laboratory 


Addressee: 
Code 1001 
CODE 1220 
Code 2628 
Code 4000 
Code 4700 


Attn: Name/Code 
T. Coffey 
Security 

TID Distribution 
W. Ellis 
S. Ossakow 


1 copy 
1 copy 


1 copy 
1 copy 


1 copy 


1 copy 
1 copy 
22 copies 
1 copy 
26 copies 



















Code 4700.1 - A.W. Ali 1 copy 

Code 4710 - C. Kapetanakos 1 copy 

Code 4740 - Branch Office 25 copies 

Code 4740 - W. Black 1 copy 

Code 4740 - A. Fliflet 1 copy 

Code 4740 - S. Gold 1 copy 

Code 4740 - A. Kinkead 1 copy 

Code 4740 - W.M. Manheimer 1 copy 

Code 4740 - M.E. Read 1 copy 

Code 4740 - M. Rhinewine 1 copy 

Code 4770 - G. Cooperstein 1 copy 

Code 4790 - B. Hui 1 copy 

Code 4790 - C.M. Hui 1 copy 

Code 4790 - Y.Y. Lau 1 copy 

Code 4790 - P. Sprangle 1 copy 

Code 5700 - L.A. Cosby 1 copy 

Code 6840 - S.Y. Ahn 1 copy 

Code 6840 - A. Ganguly 1 copy 

Code 6840 - R.K. Parker 1 copy 

Code 6840 - N.R. Vanderplaats 1 copy 

Code 6850 - L.R. Whicker 1 copy 

Code 6875 - R. Wagner 1 copy 


Naval Sea Systems Command 
Department of the Navy 
Washington, D.C. 20362 
Attn: Commander George Bates 

PMS 405-300 1 copy 

Northrop Corporation 
Defense Systems Division 
600 Hick3 Rd. 

Rolling Meadows, Illinois 60008 

Attn: Dr. Gunter Dohler 1 copy 









Optical Sciences Center 
University of Arizona 
Tucson, Arizona 85721 

Attn: Dr. Willis E. Lamb, Jr. 1 copy 

OSD/SDIO 

Attn: 1ST (Dr. H. Brandt) 

Washington, D.C. 20301-7100 1 copy 

Pacific Missile Test Center 
Code 0141-5 

Point Muga, California 93042 

Attn: Will E. Chandler 1 copy 

Physical Dynamics, Inc. 

P.0. Box 10367 

Oakland, California 94610 

Attn: A. Thomson 1 copy 

Physics International 

2700 Merced Street 

San Leandro, California 94577 

Attn: Dr. J. Benford 1 copy 

Princeton Plasma 
Plasma Physics Laboratory 
James Forrestal Campus 
P.0. Box 451 

Princeton, New Jersey 08544 
Attn: Dr. H. Hsuan 
Dr. J. Doane 

Quantum Institute 
University of California 
Santa Barbara, California 93106 
Attn: Dr. L. Elias 1 copy 

Raytheon Company 

Microwave Power Tube Division 

Foundry Avenue 

Waltham, Massachusetts 02154 

Attn: N. Dionne 1 copy 

Sandia National Laboratories 
0RG. 1231, P.0. Box 5800 
Albuquerque, New Mexico 87185 


Attn: Dr. Thomas P. Wright 1 copy 

Mr. J.E. Powell 1 copy 

Dr. J. Hoffman 1 copy 

Dr. W.P. Ballard 1 copy 

Dr. C. Clark 1 copy 


2 copies 
1 copy 











Science Applications, Inc. 

1710 Goodridge Dr. 

McLean, Virginia 22102 

Attn: Adam Drobot 1 copy 

P. Vitrello 1 copy 

Stanford University 

High Energy Physics Laboratory 

Stanford, California 94305 

Attn: Dr. T.I. Smith 1 copy 

TRW, Inc. 

Space and Technology Group 
Suite 2600 

1000 Wilson Boulevard 
Arlington, VA 22209 

Attn: Dr. Neil C. Schoen 1 copy 

TRW, Inc. 

Redondo Beach, California 90278 
Attn: Dr. H. Boehmer 
Dr. T. Romisser 

University of California 
Physics Department 
Irvine, California 92717 
Attn: Dr. G. Benford 
Dr. N. Rostoker 

University of California 
Department of Physics 
Los Angeles, CA 90024 
Attn: Dr. A.T. Lin 

Dr. N. Luhmann 
Dr. D. McDermott 

University of Maryland 
Department of Electrical Engineering 
College Park, Maryland 20742 
Attn: Dr. V. L. Granatstein 

Dr. W. W. Destler 

University of Maryland 
Laboratory for Plasma and Fusion 


Energy Studies 

College Park, Maryland 20742 

Attn: Dr. Jhan Varyan Heilman 1 copy 

Dr. John McAdoo 1 copy 

Dr. John Finn 1 copy 

Dr. Baruch Levush 1 copy 

Dr. Tom Antonsen 1 copy 

Dr. Edward Ott 1 copy 


1 copy 
1 copy 


1 copy 
1 copy 
1 copy 


1 copy 
1 copy 


1 copy 
1 copy 










University of Tennessee 
Dept, of Electrical Engr. 

Knoxville, Tennessee 37916 

Attn: Dr. I. Alexeff 1 copy 

University of New Mexico 
Department of Physics and Astronomy 
800 Yale Blvd, N.E. 

Albuquerque, New Mexico 87131 

Attn: Dr. Gerald T. Moore 1 copy 

University of Utah 

Deparment of Electrical Engineering 

3053 Merrill Engineering Bldg. 

Salt Lake City, Utah 84112 

Attn: Dr. Larry Barnett 1 copy 

Dr. J. Mark Baird 1 copy 


U. S. Naval Academy 
Annapolis, Maryland 21402-5021 


U. S. Army 
Harry Diamond Labs 
2800 Powder Mill Road 
Adelphi, Maryland 20783*1145 
Attn: Dr. Edward Brown 

Dr. Michael Chaffey 


Varian Associates 
611 Hansen Way 
Palo Alto, California 


94303 


Attn: 


Dr. 

Dr. 

Dr. 

Dr. 


H. Jory 
David Stone 
Kevin Felch 
A. Salop 


Yale University 
Applied Physics 
Madison Lab 
P.0. Box 2159 
Yale Station 

New Haven, Connecticut 06520 
Attn: Dr. N. Ebrahim 

Dr. I. Bernstein 


Director of Research 
L’.S. Naval Academy 
Annapolis, MD 21402 


2Copies 


Code 1220 


1 Conv 


31 


1 copy 


1 copy 
1 copy 


1 copy 
1 copy 
1 copy 
1 copy 


1 copy 
1 copy 














£N V 

f ILmeD 

II 

DT It- 





